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We show here that orbital degrees of freedom produce a distinct signature in the magnetic exci- 
tation spectrum of iron-based superconductors above the magnetic ordering temperature. Because 
d xz and d yz orbitals are strongly connected with Fermi surface topology, the nature of magnetic 
excitations can be modified significantly due to the presence of either static or fluctuating orbital 
correlations. Within a five-orbital itinerant model, we show that static orbital order generally leads 
to an enhancement of commensurate magnetic excitations even when the original Fermi surface lacks 
nesting at (n,0) or (0, n). When long-range orbital order is absent, Gaussian fluctuations beyond 
the standard random-phase approximation capture the effects of fluctuating orbital correlations on 
the magnetic excitations. We find that commensurate magnetic excitations can also be enhanced if 
the orbital correlations are strong. Our results offer a natural explanation for the incommensurate- 
to-commensurate transformation observed in a recent neutron scattering measurement (Z. Xu, et. 
al., arXiv:1201.4404 ), and we propose that this unusual transformation is an important signature to 
distinguish orbital from spin physics in the normal state of pnictides. Implications for the magnetic 
and superconducting states are discussed. 



PACS numbers: 

I. INTRODUCTION 

The orbital-dependence of Fermi surface pockets due 
to the partial occupancy of multiple d-states of iron is 
a key microscopic feature which distinguishes iron-based 
superconductors from cuprates. It is natural, therefore, 
that both orbital and spin-based scenarios have been ad- 
vanced to account for the structural and magnetic transi- 
tions in these materials. In the orbital-based scenario^— , 
the structural phase transition is induced by a sponta- 
neous orbital order in the quasi-one-dimensional d xz and 
d yz orbitals which breaks the C4 symmetry. Such order 
induces a stripe-like antiferromagnetic (AFM) transition 
at a lower temperature. 

In contrast, some have proposed that a spin-only 
model accounts for both structural and magnetic tran- 
sitions in pnictides. While the source of the interac- 
tions that drive the magnetism differ in the strong- 
coupling^— and weak-couplingiiiS. models proposed, 
both spin-based scenarios eventually end up with the 
same effective theory^ with a magnetic order parame- 
ter lying in the manifold of Z2 x 0(3). The structural 
phase transition is interpreted as a breaking of the Z2 
symmetry, and the subsequent stripe-like AFM transi- 
tion breaks the 0(3) symmetry. The main difficulty in 
distinguishing between these two scenarios is that both 
lead to the same Ginzburg-Landau effective action. 

Although intensive research efforts has focused on the 
spin excitation spectrum in stripe-like AFMH~— and 
superconductingiZr— states, a theoretical study on the 
effect of orbital order and fluctuations on the spin ex- 
citation spectrum in the paramagnetic normal state is 
far from complete. It is this problem that we tackle 



here. Our line of reasoning is based on the fact that 
if the spin interaction is the only energy scale driving 
both phase transitions, the difference in magnetic exci- 
tations between the normal and ordered states should be 
tied solely to a redistribution of spectral weight. The 
positions of the peaks in the Brillouin zone, namely the 
magnetic spectrum, should not undergo any qualitative 
change. On the other hand, if there is an additional en- 
ergy scale playing an equally important role, then this 
extra energy scale could result in qualitatively distinct 
changes in the magnetic spectrum above and below the 
ordering transition. 

In this paper, we analyze the magnetic excitations 
within a five-orbital itinerant model with generalized 
on-site Hubbard interactions, and superconducting^— 
states, Since our main purpose is to identify several fea- 
tures which can not be captured without involving or- 
bital degrees of freedom, an orbital-dependent non-spin 
interaction is introduced in order to investigate system- 
atically how orbital correlations affect the magnetic exci- 
tations in the paramagnetic normal state. We find that 
orbital order always enhances one of the commensurate 
spin excitation wave vectors even when the original Fermi 
surfaces lack nesting at these commensurate wave vec- 
tors. Moreover, to include the effect of fluctuating or- 
bital correlations on the magnetic excitations, we derive 
a Gaussian fluctuation model to go beyond the standard 
random-phase approximation (RPA) treatment. Even in 
this case, we find that fluctuating orbital correlations also 
tend to shift the spectral weight toward the commensu- 
rate wave vector, resulting in a broad spectrum near the 
commensurate wave vectors as seen in neutron scattering 
measurements. We show that our results offer a natural 
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explanation for the incommensurate-to-commensurate 
transformation observed in the inelastic neutron scat- 
tering measurement on Fei-xNi^Teo.sSeo.g^, which can 
hardly be explained by spin-based scenarios. As there is 
no reason to suspect that one model is insufficient to ex- 
plain the normal state of the pnictides, we conclude that 
orbital physics is the key player in driving both the struc- 
tural as well as the magnetic transitions in the iron-based 
superconductors. We propose further experimental stud- 
ies on magnetic excitations at higher temperature in the 
normal state to further pin down the orbital scenario. 



II. FORMALISM 

A. Random phase approximation (RPA) theory for 
magnetic excitation 

Starting from a tight-binding model in the unfolded 
Brillouin zone with one Fe atom per unit cell, we write 
the model Hamiltonian as 



H — H t + Hqo + Hi 



Hi = ^2Un iat n ial +^T(U'- ^)n ia n ib 



i,b>a 



2 JS ia ■ S. lb + J' (p ia p\ b + h.c.) , (1) 

i,b>a 



where 
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C — J - 

Via — C-i a ^Ci a ^. 
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We adopt the relationships U' = U — 2 J, J = J' and 
throughout the paper we use the values of U — 2.0eV 
and J = 0.2eV. H t is the five-orbital tight-binding model 
fitted by Graser et, alJ&, and Hoo is the effective inter- 
action for 'orbital order' which we put in by hand as 



E 

k.a 



i c L,Ak + Q) c ^,a(k) - cl z>u (k + q)c yz>a {k)) . 



(3) 



We do not justify the microscopic origin of Hqo since the 
conclusion is independent of the mechanism for orbital 
order. However, we will comment on the differences in 
the magnetic excitations for different scenarios in the dis- 
cussion section. For ferro-orbital order, we have rj q = —r/a 
for q = 0, which leads to a mean-field Hamiltonian of the 
form, 



H MF = H t 



{ c xz,a(k) c xz,a-(k) C yz,cr {^) C yz,a (&)) 



fe . a 



with 4> = ^(C<f=o) solved self-consistently. 

It is convenient to use the basis, ipt = 

(ct , ct , ct , ct . ct ) , and 

v k,xz,cr k,yz,cr k,xy,a k,x z —y z ,a k,3z 2 —r 2 ,a ' 

introduce U? such that (U?) 1 H MF (k,a)U? = 

diag. [Et x , • • • , Et g ] . It then follows that the 
bare response function can be expressed as 



b.m 



(5) 



Using the same notation as Graser et. aZ.— , we express 
the spin response function to one-loop as 

X RPA (q, W, 4>) = X (0) (q, W, 4>) [l - V S X (0) («f, w, 4>)} ~* (6) 
where V s is the spin interaction kernel derived in Ref. [l9[ . 



B. Gaussian fluctuations in the absence of 
long-range order 

Even above the transition for long-ranged orbital or- 
der, the magnetic excitations can still be influenced by 
the fluctuations in the proximity to the orbitally ordered 
state. In this section, the minimal Gaussian fluctua- 
tion model is derived from the path integral formalism. 
Since our focus is on the interplay between spin-flip exci- 
tations and orbital ordering in the normal state with- 
out any long-range order, we introduce the Hubbard- 
Stratonovich fields to decouple Hqo an d the spin-flip 
channels in Hi 



(7) 



and drop all the other terms. The partition function 
takes the form, 

Z = J D[c}D[J}e- So - Sl - So ° 

D[(j)]D[(j)*}D[M]D[M*]D[c]D[c^]e~ So ~ SM+s ' 00 



Sm = ^2 ^2[M ab (q,iuj n )S cd (q,iuj n ) 

%,<T,iu n abcd 

+ M ab (-q,-iuj n )S^ d (q,iuj n ) 
+ \y S ]ab,cd(& iun)M ab (-q, -iu n )M cd (q, iw n )] 
S oo = E [<K9> « w nR(-9, -iuj n ) + <j}(-q, -iu n )((q,iuj n ) 



q,lU) n 
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Because we are interested in the case where the system 
is on the verge of ferro-orbital order but not close to 
any instability in spin channel, the largest weight in the 
partition function comes from the Gaussian fluctuations 
around (q — 0, itu n = 0) in the <j> field. Therefore, we keep 
only the <fi(q = 0, iu n = 0) field in Eq. [8] and treat the 
spin part by the standard saddle point approximation. 
After integrating out the fermionic fields and expanding 
all the Hubbard-Stratonovich fields up to quadratic or- 
der, we arrive at a partition function of the form, 



(4=0 meV 



D[<l)(0)]D[<l>(0)*]D[M}D[M*]e 

xM ab (-q, -iu n )M cd (q, iu n )] 
S oa = A- 2 0*(O)0(O), 



-S m -S oa 



(9) 



where A = — H(q = 0, iu = 0) is the stiffness for 
4>(0) field renormalized by the bubble diagrams. Then it 
is straightforward to arrive at the spin response function 
with minimal Gaussian fluctuations: 

x (?, «) jd^PJv {10 > 

where A 2 is an effective parameter measuring the strength 
of the fluctuating orbital correlations. The physical 
meaning of the above Gaussian fluctuation model is clear. 
Because the system is on the verge of ferro-orbital order, 
the fluctuations of the (j> fields at (q = 0,iu> n = 0) are 
dominant. Either quantum or thermal fluctuations can 
create a temporary ferro-orbital order with a probabil- 
ity of e -1 ^ ) / A , which is the leading fluctuating orbital 
correlations in this case. It is worth mentioning that the 
orbital order parameter ((f) = is strictly fulfilled, so 
there is no long-range order. This is analogous to the 
formalism used by Gollub et. al. to study diamagnetism 
above T c in conventional superconductors^ 2 -. 



III. RESULTS 

Our main results are summarized in Figs. [IJ[2j and [4] 

The imaginary part of the spin susceptibility is plotted 
along the direction transvers to the Qafai — (tt, 0) for 
different cases. From Figs. HJa)-(c), it can be seen clearly 
that with the increase of the orbital order parameter, the 
peaks in the magnetic excitation gradually move from 
incommensurate to commensurate wave-vectors, which 
can be understood solely in terms of the Fermi surface 
topology. Generally speaking, the realistic Fermi sur- 
faces of the iron-based superconductors does not have 
perfect nesting at (ir, 0) and (0, 7r), and previous calcula- 
tions have shown that in the normal state, the peaks in 
the spin susceptibility are usually incommensurate away 
from (it, 0), (0, 7r) 16 i 18 i 21 i 23 . However, due to the presence 
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FIG. 1: The imaginary part of the spin susceptibility are 
plotted along q — (-K,q y ) for (a) <j> = 0, (b) <f) = 30 meV, 
(c) <f> = 60 meV, and (d) Gaussian fluctuations with A 2 = 60 
meV. In this figure, the chemical potential fi is set to be in 
the tight-binding model proposed in Ref. 18], corresponding 
to an electron doping of about 5%. 



of orbital order, the Fermi surfaces are distorted along 
one of the in-plane axes, which always helps the enhance- 
ment of the magnetic excitations at (w, 0) for <f> > (or 
(0, 7r) for (f> < 0). It follows then that the stripe-like anti- 
ferromagnetism accompanying pre-existing orbital order 
should have Qafm = (^,0) for (j) > or Qafai — (0, 7r) 
for (ft < 0, which has been confirmed by previous mean- 
field 2 ^ and first principle calculations 3 . 

The magnetic excitations are also strongly modified 
without long-range orbital order but with fluctuating or- 
bital correlations. In this case, because the fluctuating 
orbital correlations induce fluctuations of the Fermi sur- 
face shape, the magnetic excitations at the incommen- 
surate wave vectors are significantly diffused toward the 
commensurate ones. This results in a broadened spec- 
trum near (it, 0), which can be seen from the results of 
our Gaussian fluctuation model shown in Fig. [Tfcl) . The 
imaginary part of the spin susceptibility are plotted for 
different doping level in Fig. [2l which shows the same 
behavior. 

It is remarkable to see that our results offer a natural 
explanation for several puzzles present in previous inelas- 
tic neutron scattering measurements. As shown in Fig. 
[31 for superconducting Fei-^Ni^Teo.sSeo.s with x = 0.04, 
the overall shape of the magnetic excitations at low en- 
ergy udergoes a distinct transformation from two incom- 
mensurate vertical columns to a broad U-shaped spec- 
trum centered at the commensurate wave vectors. As 
shown in Ref. pol ], the change occurs at T onset ~ 3T C . 
Although we can not obtain the temperature-dependence 
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FIG. 2: The imaginary part of the spin susceptibility are 

plotted for doping level about 9% for (a) <f> = 0, (b) <j> = 10 FIG. 3: Inelastic neutron scattering measurement (data from 
meV, (c) <t> = 20 meV, and (d) Gaussian fluctuations with X u et. alM) of Fei^Ni^Teo.sSeo.s with x = 0.04 along the 
A =30 meV. transverse direction to Qafm for (a) high temperature normal 

state (T > Tfi uc ) and (b) the superconducting state (T < 
T c ). An incommensurate-to-conmmensurate transformation 
is observed as the sample is cooled down. 



of A 2 from a microscopic calculation within the current 
model, it is expected that A 2 should gradually increase 
if the system gets closer to the orbitally ordered state 
as the temperature is lowered. Consequently, this un- 
usual incommensurate-to-commensurate transformation 
can be captured by our theory, as shown in Fig. [U if we 
identify T onset as the onset temperature of the fluctuating 
orbital correlations Tfi uc . Moreover, a weak lattice dis- 
tortion has been found in x-ray measurements 2 ^ 2 ^, which 
gives further support for our orbital-based explanation. 

Another interesting observation pointed out by Ewings 
et. al. 2 ~ in inelastic neutron scattering measurement for 
SrFe2 As2 is that after a thorough comparison between ex- 
perimental data and theories, the multiorbital itinerant 
model treated with RPA- captures more details of the 
experimental data compared to a local moment model, 
except the theory obtained an incommensurate spectrum 
in the normal state which is not seen in experiments. 
This inconsistency can be easily resolved by the inclu- 
sion of the Gaussian fluctuation model presented in this 
paper. As seen in Fig. 2J while the incommensurate 
peaks for T > Tfi uc similar to Ref. !16| are reproduced 
in our calculations, the spectral weights of these incom- 
mensurate peaks are gradually shifted to (n, 0) as the 
static and fluctuating orbital correlations are taken into 
account. Another subtle point worthy of mentioning is 
that there is no soft mode at q = (it, n) in our results, 
which is an important advantage of itinerant model over 
local moment model as mentioned in Ref. 1211] . 



IV. DISCUSSION 

First, we would like to discuss the nature of the mag- 
netic excitations in spin-based scenarios proposed for 
iron-based superconductors. For the local spin scenario 
with frustrating J\-Ji interactions, it is generally hard 
to obtain peaks at incommensurate wave vector in the 
spin excitation spectrum in the paramagnetic state. One 
possibility is to include more longer-range spin interac- 
tions (e.g., Jz^. Even so, it would be difficult to under- 
stand why the system would start from an incommensu- 
rate spin excitation spectrum at high temperature and 
then evolve into a stripe-like AFM with a commensu- 
rate wave vector. For the SDW scenario, there are only 
incommensurate magnetic excitations in the high tem- 
perature normal state, as seen in Ref. [2fj| . and hence 
this model is not applicable. As a result, we conclude 
that the incommensurate-to-commensurate transforma- 
tion in the magnetic excitations is a unique signature fa- 
voring orbital-based over spin-based scenarios. However, 
we also recognize the fact that the size of the instanta- 
neous moments is too large to be described by a purely 
itinerant model in all current existing neutron scattering 
data. This suggests that one needs to take into account 
the orbital correlations emerging from an itinerant model 
to obtain the spin interactions correctly. Additional in- 
elastic neutron scattering measurements in the high tem- 
perature normal state are certainly necessary to resolve 
the long-standing controversy between the orbital-based 
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FIG. 4: Demonstration of the evolution of the imaginary part 
of the spin susceptibility along q = (tt, q y ) as the tempera- 
ture is lowered from high temperature. The parameters in 
the tight-binding model are the same as Fig. [T] and a broad- 
ening in the g-space of qo = 0.17T is introduced with the form 

of Im X B (g» = Y,„ e- (? - 9 " )2/9 °Irn X (?>)/i:< ? e^"^^ 
for the ease of comparison with the experimental data in Fig. 
[3] For T > Tfi uc , orbital fluctuations do not affect the mag- 
netic excitations. In this case, Hoo is completely turned off. 
For Tfi uc > T > T s , the system does not have long-range 
orbital order but rather has fluctuating orbital correlations, 
which can be described at the Gaussian level (in the plot, 
A 2 = 30 meV is used). For T < T s , an orbital order is formed 
and the magnetic excitation is commensurate. The orbital 
order parameter in the plot is <j> — 20 meV. 



and spin-based scenarios. 

Second, we comment briefly on materials. FeSe is 
a known iron-based supercondcutor exhibiting only a 
structural transition and no magnetic transitio n ' 28 and 
hence is a promising prototype for the orbital-based sce- 
nario as well as the incommensurate-to-commensurate 
transformation in its magnetic excitation spectrum. The 
families of iron-based superconductors whose structural 
and magnetic transition temperatures are well-separated, 
including LaOFeAs^, CeOFeAs2&, NaFeAs 3 ^, etc., 
should also be good candidates to realize the physics out- 
lined in this paper. There have been neutron scattering 
measurements for parent compounds of 122, for exam- 
ple, BaFe 2 As2^, CaFe 2 As2^^, and SrFe 2 As2^-, and in- 
terestingly they show very different trends. While the 
magnetic excitations in CaFe 2 As 2 can be well-described 
by the local spin scenario both in the magnetic and the 
normal states, itinerant models seems to work better in 
BaFe 2 As 2 and SrFe 2 As 2 . As suggested in Ref. [2lJ, the 
reason why the local spin scenario fails in SrFe 2 As 2 is 
that in the high temperature normal state, the magnetic 



excitations at momentum (tt, it) are pushed to high en- 
ergy which clearly contradicts the prediction of the local 
spin scenario unless one accepts that the anisotropic spin 
interactions can still be present even above the magnetic 
transitions. The most serious problem with the itiner- 
ant model for SrFe 2 As 2 is the prediction of an incom- 
mensurate excitation spectrum in the high temperature 
normal state. This can be resolved with the Gaussian 
fluctuating orbital correlations proposed here. In other 
words, in order to explain the neutron scattering data 
of SrFe 2 As 2 , the orbital degrees of freedom are neces- 
sary even in a spin-based scenario. Therefore, since the 
previous work on BaFe 2 As 2 only measured up to 150K, 
slightly above the magnetic transition temperature (~ 
140K), it would be essential to perform additional neu- 
tron scattering measurement at much higher tempera- 
ture to see whether the features observed in SrFe 2 As 2 ap- 
pear in BaFe 2 As 2 . Should they be observed, this would 
further substantiate the orbital-based physics proposed 
here. 

Third, we discuss the implication of our theory for su- 
perconductivity. Since most iron-based superconductors 
require doping away from the parent compounds to in- 
duce superconductivity, this strongly suggests that the 
Fermi surfaces of superconducting samples generally do 
not have perfect nesting at (n,0) and (0,tt). Conse- 
quently, if the pairing mechanism of iron-based supercon- 
ductors is really through spin fluctuation o 18 ' 23 , the neces- 
sary commensurate spin excitation in this mechanism can 
not be obtained without the fluctuating orbital correla- 
tions discussed in this paper. This implies that the role of 
the orbital degrees of freedom in the superconductivity 
in iron-based superconductors might be very profound. 
Recently, orbital-dependent superconducting gaps have 
been observed in the superconducting Bai_ :E K; I ;Fe 2 As 2 
via Laser ARPES measurement 3 ^. The electronic struc- 
ture of the vortex core in the FeSe superconductors shows 
a strong anisotropy 3 ^, which has been interpreted as 
a consequence of competition between nematicity (or- 
bital order) and superconductivit y 40 ! 41 . A very recent 
ARPES measurement 4 ^ revealed that the doping depen- 
dence of the band renormalizations of d xz and d yz or- 
bitals is strongly correlated with the enhancement of 
the superconductivity. These findings strongly suggest 
that only an inceptive understanding of superconductiv- 
ity would be obtained without including the orbital de- 
gree of freedom. We propose that the incommensurate- 
to-commensurate transformation in the magnetic excita- 
tions should be a general feature in most superconducting 
iron pnictides. 

Finally, we would like to discuss the magnetic exci- 
tation along the direction longitudinal to Qafm- Our 
model showed that the magnetic excitation still disperses 
along this direction, despite the fact that the spectral 
weight could be different from that along the transverse 
direction. Experimentally, it has been found that the 
magnetic dispersion is absent along the longitudinal di- 
rection for FeTe/Se system o 20 i 43 , which remains a puzzle 
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and is beyond our current model. One possibility is to 
exploit an effective degenerate double exchange model 5 
whose spin interactions are orbital dependent so that the 
effects of the orbital correlations discussed in this paper 
could be correctly incorporated. This work is currently 
under development. 

V. CONCLUSION 

In this paper, we have studied magnetic excitations 
in the normal state of iron based superconductors. Be- 
cause the orbitals of d xz and d yz are closely related to 
the shape of the Fermi surface, the orbital order, which 
distorts the Fermi surfaces, can have a profound effect 
on the magnetic excitations. We do not specify the ori- 
gin of the orbital ordering and fluctuations since the re- 
sults are generally valid regardless of the origin. Fur- 
thermore, we assume that the ferro-orbital fluctuations 
are the dominating ones as the structure phase transi- 
tion is approached from higher temperature. We have 
shown that orbital order always enhances one of the 
commensurate spin excitation wave vectors even when 
the original Fermi surfaces do not have a nesting at 
these commensurate wave vectors. Furthermore, even 
when long-range orbital order is absent, Gaussian fluc- 



tuations of the orbital order can still couple to the mag- 
netic excitations. We found that these Gaussian fluc- 
tuations also tend to diffuse the spectral weights to- 
ward the commensurate wave vector, resulting in a broad 
spectrum near the commensurate wave vectors as seen 
in neutron scattering measurements^. We have shown 
that this orbital-based scenario offers a natural explana- 
tion for the incommensurate-to-commensurate transfor- 
mation observed in Fei-^Ni^Teo.sSeo.s, which can hardly 
be explained by spin-based scenarios. We propose that 
this unusual transformation is an important signature 
that can distinguish orbital from spin-based physics in 
pnictides. Our theory can be tested by additional neu- 
tron scattering experiments at higher temperatures. 
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